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We highlight the remarkable evolution in the cosmic microwave background 
(CMB) power spectrum Ci as a function of multipole I over the past few years, and 
in the cosmological parameters for minimal inflation models derived from it: from 
anisotropy results before 2000; in 2000, 2001, from Boomerang, Maxima and the De- 
gree Angular Scale Interferometer (DASI), extending £ to approximately 1000; and 
in 2002, from the Cosmic Background Imager (CBI), the Very Small Array (VSA), 
ARCHEOPS and the Arcminute Cosmology Bolometer Array Receiver (ACBAR), 
extending I to approximately 3000, with more from Boomerang and DASI as well. 
Pre-WMAP (pre- Wilkinson Microwave Anisotropy Probe) optimal bandpowers are 
in good agreement with each other and with the exquisite one-year WMAP re- 
sults unveiled in February 2003, which now dominate the I < 600 bands. These 
CMB experiments significantly increased the case for accelerated expansion in the 
early universe (the inflationary paradigm) and at the current epoch (dark energy 
dominance) when they were combined with 'prior' probabilities on the parameters. 

The minimal inflation parameter set, {ui,, Lo c d m , Fltot, n s , Tc, os}, is applied 
in the same way to the evolving data. Ci database and Monte Carlo Markov Chain 
methods are shown to give similar values, highly stable over time and for different 
prior choices, with the increasing precision best characterized by decreasing errors 
on uncorrelated 'parameter eigenmodes'. Priors applied range from weak ones to 
stronger constraints from the expansion rate (HST-h), from cosmic acceleration 
from supernovae (SN1) and from galaxy clustering, gravitational lcnsing and local 
cluster abundance (LSS). After marginalizing over the other cosmic and experi- 
mental variables for the weak+LSS prior, the pre-WMAP data of January 2003 
compared with the post-WMAP data of March 2003 give fi tot = 1.03±o °4 com- 
pared with 1.02t°;^, consistent with (non-baroque) inflation theory. Adding the 
flat Qtot — 1 prior, we find a nearly scale invariant spectrum, n s = 0.95^q'q4 com- 
pared with 0.971qq2- The evidence for a logarithmic variation of the spectral tilt 
is < 2cr. The densities are: for baryons, uit = fif,h 2 = 0.0217+q'qo2 compared with 
0.0228to;°°i, near the the Big Bang nucleosynthesis estimate of 0.0214 ± 0.002; for 
CDM, uj cdm = n cdm h 2 = 0.126toml compared with 0.121±° ^°; and for the sub- 
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stantial dark (unclustered) energy, Q\ 0.66l ; 09 compared with 0.70i 05 . The 
dark energy pressure-to-density ratio wq is not well constrained by our weak+LSS 
prior, but adding SN1 gives wq < —0.7 for January 2003 and March 2003, consis- 
tent with the wq=—1 cosmological constant case. We find as — 0.89^^07 compared 
with 0.86to'o4> implying a sizable Sunyaev-Zeldovich (SZ) effect from clusters and 
groups; the high £ power found in the January 2003 data suggest erg ~ 0.94+0.16 is 
needed to be SZ-compatible. 

Keywords: cosmology, particle physics, dark matter, dark energy 

1. The evolution of CMB spectra and cosmic parameters 

We have been in the midst of a remarkable outpouring of results from the CMB since 
1999. The Royal Society Discussion Meeting focused on the eight pre- Wilkinson Mi- 
crowave Anisotropy Probe (pre-WMAP) announcements made in 2002. The WMAP 
release was three weeks later, and a before- WMAP discussion without an after- 
WMAP discussion is unthinkable now. This paper applies the same methods of 
analysis to WMAP as to the earlier CMB experiments to put its singular forward 
step into context. In § 2, we describe the different experiments that have contributed 
to the evolving picture. For more background material on methods and references, 
see Bond (1996), Bond & Crittenden (2001), Bond et al. (2003a,b), Sievers et al. 
(2003), Ruhl et al. (2003) and Goldstein et al. (2003). 

(a) Grand unified spectra compared with WMAP 

Figure 1 shows how the CMB power spectrum, C e = 1(1 + l)(\T tm \ 2 ) / (2ir) de- 
fined in terms of CMB temperature anisotropy multipoles Ti m , changed from the 
pre-WMAP determination shown at the Royal Society meeting using the data to 
January 2003 to the post- WMAP determination of March 2003; the agreement is 
excellent. How we got there is shown in figure 2, providing snapshots compared with 
WMAP for data that was available in January 2000, January 2002, June 2002, as 
well as January 2003. Accompanying this story is a convergence with decreasing 
errors over time on the values of the cosmological parameters given in figures 3, 4 
and table 4. 

Optimal spectra and their error matrices are calculated in exactly the same way 
that cosmological parameters are, with the parameters now the bandpowers Cb in 
the chosen €-bins, b. Additionally characterizing each experiment there are calibra- 
tion uncertainties and often beam uncertainties, each adding additional parameters. 
Sample values for these are given in § 3b. 

(b) Basic parameters characterizing the early Universe and CMB transport 

Our philosophy has been to consider minimal models first, then see how progres- 
sive relaxation of the constraints on the inflation models, at the expense of increas- 
ing baroqueness, causes the parameter errors to open up. We adopt the basic set 
of seven cosmological parameters {u>b, w c rf m , fife, n s , tc, ln'P<i>(fc n )} to facilitate 
comparison with results in Lange et al. (2001), Jaffe et al. (2001), Netterfield et 
al. (2002), Sievers et al. (2003) and Goldstein et al. (2003). How the values have 
converged upon the bull's-eye 2a determinations with WMAP is shown in figure 3. 
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In spite of the great success in extending the spectrum to high £, the evolution of 
the parameter errors was not that strong after January 2002 until WMAP. This 
is because the Ci model space is restrictive for inflation-based models, with high 
I intimately related to lower £. On the other hand, when the experiments were 
treated individually (always with COBE-DMR), their 2a contours were all circling 
the bull's-eyes (Sievers et al. 2003). 

The transport of the radiation through the era of photon decoupling is sensitive 
to the physical density of all of the species of particles present then, Wj = IX, h 2 . We 
use two parameters, u>b for baryons and io c dm for cold dark matter, to characterize 
this, but we should add u>hdm for hot dark matter (massive but light neutrinos), 
and uj er for the relativistic particles present at that time (photons, very light neu- 
trinos, and possibly weakly interacting products of late time particle decays). Here 
the latter is fixed for the conventional three species of relativistic neutrinos plus 
photons. The total matter density is 

U m = LUb + LOcdrn + <^hdm ■ 

Another two parameters characterize the transport from decoupling to the present, 
the vacuum or dark energy, encoded in a cosmological constant and the cur- 
vature energy Qfc = 1 — VL to t- Of course fifc also determines the mean geometry. 
(When one wishes to focus on what the CMB can tell us about the nature of the 
dark energy, another parameter is often added, wq = pq/pq, where pq and pq 
are the pressure and density of the dark energy. If the vacuum or dark energy is 
reinterpreted as Qq, the energy in a scalar field Q which dominates at late times, 
it would be likely to have complex dynamics associated with it. In that case, Q and 
wq would have spatial and temporal variations (except if wq = — 1, the cosmolog- 
ical constant case). Spatial fluctuations of Q are expected to leave a direct imprint 
on the CMB for small I. This complication is typically ignored, but should not be. 
It does depend in detail upon the specific model for Q.) 

In this parameter space, h = Q^jWj) 1 / 2 and the age of the Universe to are 
derived functions of the u>j, fik,A and wq. 

Another parameter is the Compton 'optical depth' tc from a reionization red- 
shift z re h to the present, 

T C ~ 0.12(w b /0.02)(w m /0.15)- 1 / 2 ((l + z roh )/15) 3 / 2 . 

As long as tc is not too large, Ce is suppressed by a factor exp[— 2rc] on scales 
smaller than the horizon at z ro h- For typical models of hierarchical structure for- 
mation, we expect tc ^ 0.3. At the moment, even with WMAP, the CMB total 
anisotropy (TT) alone does not give such a constraint. It is the cross correlation 
of total anisotropy with polarization (TE) that leads to a detection (Kogut et al. 
2003). 

Two parameters characterize the early universe primordial power spectrum of 
gravitational potential fluctuations one giving the overall power spectrum am- 
plitude "P$(fc„), and one defining the shape, a spectral tilt 

n s (k n ) = 1 + d In V^/d In k , 

at some (comoving) normalization wavenumber k n . Instead of lnV®(k n ), which is 
appropriate for connecting to early universe physics, we use as a basic amplitude 
variable InCio when connecting to CMB, and In of when connecting to LSS. 
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To characterize inflation, even in the simplest models, we really need at least 
another two parameters, Vawikn) and n t (fc„), associated with the gravitational 
wave (GW) component. In inflation, the amplitude ratio Vow /T 3 * is related to n t 
to lowest order, with 0{n s —n t ) corrections at higher order (e.g. Bond 1996). There 
are also useful limiting cases for the n s — n t relation. With January 2003 data, and 
even with WMAP, the data are not powerful enough to determine much about the 
GW contribution, e.g. the WMAP team estimate the gravitational wave (tensor) 
contribution to be less than 0.72 of the scalar component in amplitude at the 95% 
CL. 

As one allows the baroqueness of the inflation models to increase, one can en- 
tertain essentially any power spectrum. This implies a fully fc-dependent n s (k) if 
one is artful enough in designing inflaton potential surfaces. The simple model 

n s (k) = n s (k n ) + [dn s {k n ) / d\nk] ln(fc/fc„) 

adds a logarithmic running index about a pivot scale k n . As figure 1 and § 4 indicate, 
this improves the fit to the data. It is also expected in inflation models, it is just a 
question of the size of the correction. 

The tensor index n t {k) could also be a function, although it does not have as 
much freedom as n s (k) in inflation. For example, it is difficult to get nt{k) to be 
positive. One can also have more types of modes present, e.g. scalar isocurvature 
modes. The data have shown for quite a while that these would have to be sub- 
dominant relative to the scalar curvature modes, and would have to be even more 
so now. 

Each experiment also contributes a parameter describing the uncertainty in the 
calibration, and possibly another for the uncertainty in the beam size. 

(c) CMB analysis pipelines: bandpowers to cosmic parameters 

In Gaussian models defined by a parameter set {y a }, the probability distribution 
of the primary anisotropics is fully encoded in the isotropic power spectrum Ce(y a ) 
- as long as there is no preferred orientation (as might occur for small universes that 
are topologically nontrivial). The observed bandpowers for an individual experiment 
can then be tested against theoretical bandpowers Cb(y a ), which are averages of the 
C^s over ^-space 'window functions' ipbe appropriate to the bands for the experiment 
in question. This represents a huge compression of the entire dataset and makes 
large model space computations feasible. 

To use this information to estimate cosmological parameters, the entire likeli- 
hood surface as a function of the {Cb} is needed with sufficient accuracy that the 
parameter estimations are not biased. It has been shown that individual bandpowers 
Cb have distributions well characterized by a lognormal distribution in the variable 
Cb + Cjvb, where Cjvb is an estimate of the noise in the band (Bond, Jaffe & Knox 
2000). The coupling between bandpowers is included as a weak correction, relying 
on the band-to-band correlations being relatively small — a demand imposed in 
the data analysis phase. What comes out are entropies, S(y a ), i.e. log likelihoods. 
A slightly modified version of this prescription is used for WMAP (Verde et al. 
(2003); see also § 2e.) 

There are two approaches to sampling the set {y a } that we have used. The 
main workhorse throughout our analyses up to January 2003 used fixed grids: a 
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discrete set of parameter values are chosen a priori for six of the seven cosmic 
variables, with spacings in each of the dimensions designed by hand to be adaptively 
concentrated about the most probable values, but with sufficient spread to ensure 
that tails and multiple solution regions are well explored. The current database for 
the 'minimal inflation' parameter model contains 8.5 x 10 6 models, with dimensions 
f 5x f3x f 5x 12 x 3f x f 1 for the 'external parameter' set {uJb, uj c d m , ^tot, ^A,n S7 rc}, 
with edge cutouts requiring Q TO > 0.1. 

The seventh (amplitude) parameter, InCio or hicf, and the experimental cali- 
bration and beam uncertainty variables are continuous. They relax to their maximum- 
likelihood values, with errors characterized by the second derivative of the likelihood 
function. The number of these 'internal' continuous parameters may become much 
larger if we split the amplitude parameter into many, one for each band in ^-space 
(or fc-space if three-dimensional power spectra are the target). The shape of 
an assumed spectrum multiplies the adopted window functions for the bands. For 

(s) 

the optimal bandpowers that combine experiments together in figure 2, C\ is usu- 
ally varied to test robustness of the results, but an ensemble of external parameter 
models can be applied, e.g. in broken scale-invariance applications in /c-space. 

The first stage output is large entropy files that include maximum-likelihood 
values and Fisher matrices for the internal variables. These files are then picked 
up in postprocessing as various prior probabilities are applied, marginalizations are 
done, and one-dimensional (ID) and two-dimensional (2D) statistics are computed. 

An advantage of a fixed grid is that it has allowed us to quickly check many 
prior cases and many experimental combinations, all on the same footing. Calibra- 
tion uncertainties are handled cither at the entropy stage with the complete experi- 
mental mix, or in postprocessing, since we know the amplitude distributions. When 
analysing an experiment, these operations are done again and again, as different 
hypotheses, band widths and positionings, estimation techniques, source removal 
methods, etc., are applied to the bandpowers for the experiments in question, so 
speed is essential. Combining DMR with the experiments was always a first step; 
now WMAP takes that role. 

As new experiments are added which are qualitative improvements (like WMAP), 
errors may become smaller than grid spacings and further adaptivity of the grid is 
needed. This is so even with good interpolations and smoothing. We set a floor on 
parameter errors to be half the grid spacing of the encompassed grid: a value that 
was never reached before WMAP but has been with WMAP for a few parameters 
strongly bundled into the top few parameter eigenmodes and some priors. 

The second approach is the Monte Carlo Markov Chain (MCMC) method 
(Metropolis et al. 1953; Christensen & Meyer 2001; Christensen et al. 2001; Lewis 
& Bridle 2002; Verde et al. 2003). It develops a set of independent chains, each a 
small (unstructured) grid on the parameter space that is constructed 'on the fly' 
rather than a priori. The elements of the chains are sampled according to well- 
developed MCMC algorithms designed to make the next step independent of prior 
ones. The spacing of the models computed changes with experimental combinations 
and priors adopted. As with the fixed-grid methods, some priors can be applied in 
postprocessing, which speeds up the procedure. These MCMC methods have now 
become feasible for CMB analyses because the Ce(y a ) computations with CMB fast 
(Seljak & Zaldarriaga 1996) or CAMB (Lewis, Challinor & Lasenby 2000) are effi- 
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cient algorithmically, and are rapid enough on large numbers of models because of 
the remarkable speedup of individual computer processors in recent years. A nice 
Fortran 90 package is publicly available to do this (Lewis & Bridle 2002). To do 
the statistics well, one needs not just many elements in a chain, but a number of 
chains. Thus it was really the advent of massively parallel machines that is allow- 
ing MCMC to become a major working tool for repetitive CMB and LSS analyses. 
For example, it is the method adopted by the WMAP team (Verde et al. 2003; 
Spergel et al. 2003), who applied it to the minimal six-parameter flat model, with 
the amplitude as an external parameter, and seven-parameter models with the Qk, 
uuq, dn s /dhik or nt allowed to vary in turn. They ran four chains and about 30000 
models per chain to define their distributions. We have adapted Cosmomc to the 
parameter choices and ranges of our Ce database, and the many prior cases and ex- 
perimental combinations used, to facilitate comparisons; e.g. tc can go out as far as 
0.7, which ensures likelihood drop-off to zero, but places sampling challenges. The 
major challenge for MCMC is to sample well the curved likelihood ridges at rea- 
sonable computational cost. Certain nonlinear combinations of our basic variables 
can help to straighten out the likelihood surfaces, in particular those with highly 
asymmetric errors, allow for more efficient and accurate computation, whether they 
be used in MCMC (Kosowsky et al. 2002; Chu, Kaplinghat & Knox 2003; Verde 
et al. 2003) or for fixed or adaptive grids. For MCMC, another approach is to use 
variance matrices from small runs to make the steps efficient in the parameter space 
(Lewis & Bridle 2002). For each experimental mix and prior, we use 16 chains run 
until convergence tests are satisfied for all of the variables. In spite of processor 
speed, the computations remain a challenge if many cases need to be run. 



{d) Weak, HST-h, SN1, LSS & t c priors 

The parameter grids are chosen to be wide relative to conceivable cosmological 
models, yet are concentrated in the maximum-likelihood regions. The MCMC chains 
are allowed to vary over wide domains, and they automatically concentrate well. An 
important issue is the prior measure we impose upon the parameter space. Implicit 
in the adoption of a given variable set is that a uniform prior probability is chosen 
in each of the variables. If a variable is not well determined this can have a big 
influence (Lange et al. 2001). 

We usually present the cosmological conclusions we draw from our analyses of 
the various CMB experiments using noncontroversial priors, ones that almost all 
cosmologists would agree to. Thus our standard weak one used in Lange et al. 
(2001) and subsequent works requires only 0.45 < h < 0.90 and to > lOGyr. The 
addition of the flat prior has also become benign, thanks to the sharpness of the 
Clk i=a determination with the CMB rather than to the predilections of inflation 
theorists. (Although a major reduction in number of database models occurs when 
the flat flk — prior is applied, it is usually applied in the postprocessing phase.) 

Data from sources other than the CMB can be incorporated as 'prior probabil- 
ities'. A stronger prior on the Hubble parameter, HST-h, uses an ft, = 0.72 ± 0.08 
Gaussian distribution (Frecdman et al. 2001). SN1 data imposes a prior in 
n k -w Q space (Riess et al.1998, Perlmutter et al. 1999a,b, 2003). The CMB data 
apparently determine to higher accuracy than light-element-abundance obser- 
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vations coupled to Big Bang Nucleosynthesis theory (Kirkman et al. 2003), hence 
applying a BBN prior is not of much interest. 

The LSS prior we use (Bond & Jaffe 1999; Lange et al. 2001; Bond et al. 
2003b) also depends upon our parameter set. An important combination is the 
wavenumber of the horizon when the energy density in relativistic particles equals 
the energy density in nonrelativistic particles: ~ 5T~ q x h^ 1 Mpc, where Y eq — 

Q m h (1.68^/wer) 1 ^ 2 - We represent the (linear) density power spectrum by a single 
shape parameter: 

T = T eq exp[-(fi 6 (l + «~ 1 (2h) 1 ' /2 - 0.06))] 

works reasonably well, to about 3% over the region most relevant to LSS; replacing T 
by r e ff = T+ (n s — 1) /2 takes into account the main effect of spectral tilt over the LSS 
wavenumber band (Bond 1996). For low redshift clusters, the abundances determine 
a combination that is roughly og^m 56 (with the degeneracy among the combination 
broken with high redshift cluster information). Weak lensing determines a similar 
combination. 

With the wealth of data emerging from the Sloan Digital Sky Survey (SDSS) 
and the 2dF Redshift Survey (2dFRS), shape is a very powerful probe. However, 
the biasing of the galaxy distribution with respect to the mass becomes an issue if 
it is scale dependent, inviting caution - and for our purposes a weakened prior over 
what the data formally show. In the future, weak lensing should allow shape and 
amplitude to be simultaneously constrained without biasing uncertainties. 

To be explicit, our prior for In of is of form <78f2^ 56 =0.47l 002 'lo 08> distributed 
as a Gaussian (first error) smeared by a uniform (tophat) distribution (second er- 
ror) . This straddles most of the values determined from weak lensing and many of 
those estimated from cluster abundances (shown in figure 5). We have also used a 
prior shifted downward by 15% to accommodate the lower values quoted for clusters 
in the literature. Our prior for the shape parameter is r c ff =0.211q'q3'1q'q8, which 
encompasses the recent SDSS and 2dFRS results as well as results from the Auto- 
mated Plate Measuring (APM) angular survey and earlier redshift surveys. Fully 
embracing the 2dFRS galaxy power spectrum with a linear bias model to relate 
it to the total density power gives a stronger shape constraint, r c ff=0.21^Q'g3'^g. 
Although some SDSS estimates are quite consistent, e.g. T c ff =0.19^q^4'^q, differ- 
ent analysis methods and different datasets give wider ranges and the estimates do 
not incorporate possible complexities in the bias model. Thus, we have adopted a 
weak-LSS as opposed to a strong-LSS prior. Since T oc uj m /h, with the improved 
CMB estimations of uo m that arose in the January 2002 (and later) datasets, the 
shape constraint now has some similarity to an h prior (§ 4). 

Values of T eq , T and as as estimated from the CMB data are given in table 4. 
They are basically compatible with the LSS priors. One can get T e ff from the L 
and n s results in the table. 

One of the most exciting results from WMAP was the TE cross-correlation of the 
E-mode of polarization and the total intensity (T) at low £, interpreted as evidence 
for a tc — 0.16 ±0.04 detection, determined with a 'model independent' method by 
Kogut et al. (2003). The detection is not nearly as strong when ensemble-averaged 
over model space for the weak prior, as described in § 4. The TE result is explicitly 
included in the Cosmomc treatment, but only TT is included in the Ci database 
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used here. To incorporate the detection, we have constructed a tc prior, chosen to 
be broader than a 0.16 ± 0.04 Gaussian, tq — ^-^X^mX^^oq m terms of Gaussian 
and top- hat errors. The MCMC results for tq we obtain when other parameters are 
marginalized is broader still, on both sides, and so we compare parameter estimates 
with and without this prior and find for most it makes little difference It does 
have an effect on marginalized amplitude determinations, in particular skewing 
somewhat the as distribution to higher values. 

Sometimes there is 'tension' between the parameters estimated from CMB-only 
and those including non-CMB priors. This is extremely important to flag, since 
poor distribution overlap leads to smaller combined errors. 



One is tempted to open up parameter space to a much larger set. There was a 
good reason for limiting the number in the pre-WMAP days: the spectra may not 
change much as the parameters vary, manifested by near-degeneracies among them. 
It is useful to disentangle the degeneracies by making linear combinations which 
diagonalize the error correlation matrix (Ay a Ay a >) , where Ay a = y a — (y a ) and the 
averages are over the probability-weighted ensemble of models. These 'parameter 
eigenmodes' (Bond 1996; Efstathiou & Bond 1999; Lange et al. 2001; Goldstein et 
al. 2003) 



The error on the eigenmode £ Q , a a , is determined by the data and the priors. (In- 
stead of AuJb and Aw C( j m , we use Awj/wj and Au) c d m /<^cdm m the combinations so 
their errors are relative and quantitatively meaningful relative to the other vari- 
ables.) 

Until the WMAP data, only four combinations of our seven could be determined 
within ±0.1 accuracy with the CMB (five with CMB+LSS), but with WMAP 
precision, for the March 2003 data, five can be determined (six with CMB+LSS), 
and two are determined to better than ±0.01. Parameter eigenmodes arising from 
the current data are discussed further in § 4b. 

Thus, WMAP precision gives us licence to open up the parameter space more. 
Here we only do this to a limited extent, by restricting ourselves to flat universes 
and replacing Qtot by wq or by dn s (k n ) / dink. 

Both MCMC and fixed-grid approaches can have difficulty when the eigenmodes 
are precisely determined. Using variables which are nonlinear combinations of the 
{y a } motivated by the eigenmodes can aid with this, e.g. one characterizing the 
peak/dip pattern (the sound-crossing scale) and one the amplitude-Tc-n s near- 
degeneracy. Both degeneracies were exploited in limiting our Ce database storage 
requirements. 

Expressing the LSS prior in terms of T + (n s — l)/2 and In crff^ 12 only (Bond et 
al. 2003b; Lange et al. 2001; Bond & Jaffe 1999) is similar in spirit to keeping only 
the best determined 'eigenmode' from the redshift surveys and from the lensing or 
cluster surveys. However, the same mechanism that gives the acoustic peaks in Ce 
leads to oscillations in the density power for large uib/uj m ; i.e. further eigenstructure 
that would be revealed with high precision shape data. Similarly extra variables such 
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as ujhdm also lead to more eigenstructure. Higher redshift observations also break 
the In crffi^ 12 near-degeneracy. 

(/) CMB pillars 

There were 'seven pillars' of the inflation paradigm that we were looking for in 
the CMB probe: 

(1) the effects of a large scale gravitational potential at low multipoles; 

(2) the pattern of acoustic peaks and dips; 

(3) damping; 

(4) Gaussianity (maximal randomness for a given power spectrum) of the pri- 
mary anisotropies; 

(5) secondary anisotropies associated with nonlinear phenomena, due to the SZ 
thermal and kinetic effects, inhomogeneous rcionization, weak lensing, etc.; 

(6) polarization, which must be there at the ca. 10% level, along with a specific 
cross-correlation with the total intensity; 

(7) anisotropies and the associated polarization induced by gravity wave quan- 
tum noise. 

At least five, and possibly six, of these have been seen. We have known about 
pillar 1 since COBE and FIRS, and found pillars 2 and 3 in the past few years, as 
discussed in § 3c, d. 

Most, but not all, inflation models predict Gaussianity of the primary CMB 
fluctuations (pillar 4). This has been demonstrated to varying degrees with COBE, 
Maxima, Boomerang, the Cosmic Background Imager (CBI) and now with WMAP 
data. All secondary anisotropies and Galactic foregrounds will be non-Gaussian, so 
care must be taken in interpreting the inevitable deviations from Gaussianity. 

The CBI, the Arcminute Cosmology Bolometer Array Receiver (ACBAR) and 
the Berkeley Illinois Maryland Array (BIMA) may have seen evidence for the ther- 
mal SZ effect, an aspect of pillar 5 (see §§ 2g and 4a). 

Polarization (pillar 6, see § 2f), has been convincingly demonstrated. First 
there was the broadband detection by DASI of EE polarization and its TE cross- 
correlation with total intensity, at levels consistent with inflation models. Then 
WMAP unveiled the TE cross-correlation spectrum to £ ~ 400. The enhancement 
at £ < 20 is the evidence for tq = 0.16 ±0.04 and an associated redshift of reioniza- 
tion z re h ~ 15. The WMAP anticorrelation in TE observed at I <~ 100 is interpreted 
as proof that the dominant component of the perturbations giving rise to this effect 
is adiabatic. 

Pillar 7 is an extreme experimental challenge, and some inflation models have 
gravity wave induced anisotropies too small for them ever be detected (see § 2f). 

2. The CMB experiments 

(a) January 2000 

We were dealing with upper limits on anisotropies until the 1990s. These limits 
were useful in ruling out broad ranges of theoretical possibilities. Then the familiar 
2 < I < 20 multipoles at the 30/iK level were revealed by COBE. This was followed 
in the years up to April 1999 by detections, and a few upper limits (ULs), at higher 
I in 19 other ground-based (gb) or balloon-borne (bb) experiments. Some predated 
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the 1992 COBE announcementin design and even in data delivery. We have the 
intermediate angle SP91 (gb), the large angle FIRS (bb), both with strong hints 
of detection before COBE, and then, post-COBE, more Tenerife (gb), MAX (bb), 
MS AM (bb), white-dish (gb, UL), argo (bb), SP94 (gb), SK93-95 (gb), Python 
(gb), BAM (bb), CAT (gb), OVRO-22 (gb), SuZIE (gb, UL), QMAP (bb), VIPER 
(gb) and Python V (gb). Most had many fewer resolution elements than the 600 
or so for COBE. One exception was SK95, which forged new ground compressing 
raw timestreams in software on to spatially extended 'pixels'. This heterogeneous 
dataset up to that time showed evidence for a peak (Bond, Jaffe & Knox 2000), 
although it was not well localized. Improved first-peak localization occurred in sum- 
mer 1999 with the Chile-based Toco experiment (Miller et al. 1999) and in Novem- 
ber 1999 with the North American balloon test-flight of Boomerang (Mauskopf 
et al. 2000). Collectively we denote the results of all of these experiments as the 
January 2000 data. It pointed to Q, tot <~ 1, but with broad errors (figure 3). 

Technical notes 

Boom-NA, Toco, QMAP, SK95, MS AM, SP94 and SP91 include quoted cali- 
bration errors; 13 bandpower detections from other experiments do not explicitly 
include them, but are often incorporated in the quoted error bars. The three band- 
power upper limits use the 'equal variance approximation' to the likelihood (Bond, 
Jaffe & Knox 2000), as opposed to the offset lognormal which is used for the rest. 
The low quadrupole is not included in the DMR bandpowers. 

(b) January 2002: Boomerang, Maxima and DASI 

In April 2000 dramatic results to I ~ 600 from Boomerang, the first CMB long 
duration balloon flight which circled Antarctica for 10.6 days in December 1998, 
were announced (de Bernardis et al. 2000; Lange et al. 2001). This was quickly fol- 
lowed in May 2000 by results from the night flight of Maxima (Hanany et al. 2000). 
Boomerang's best resolution was 10.7' ± 1.4', about 40 times better than that of 
COBE, with tens of thousands of resolution elements. (The corresponding Gaussian 
beam filtering scale in multipolc space is ~ 800.) Maxima had a similar resolution 
but covered an order of magnitude less sky. In April 2001, the Boomerang analysis 
was improved and much more of the data were included, delivering information on 
the spectrum up to i ~ 1000 (Netterfield et al. 2002). Maxima also increased its i 
range (Lee et al. 2002). 

Boomerang had six bolometers at 150 GHz and 10 other bolometers at 90, 220 
and 400 GHz. It mapped 1800 square degrees. The April 2000 analysis used only 
one channel and 440 deg 2 . The April 2001 analysis used 4 channels at 150 GHz 
and 800 deg 2 , 1.8% of the sky. These are the data used in the January 2002 and 
June 2002 sets (Netterfield et al. 2002). In December 2002, the final Boomerang-98 
analysis was given (Ruhl et al. 2003), encompassing an area 1200 deg 2 , f s k y = 2.9%, 
using 3.5' pixels and the 150 GHz channels. This is part of the January 2003 set. 

Boomerang had a successful second long-duration balloon flight in January 2003, 
taking data for 13 days. The forecasts are for a well determined TT power spectrum 
even in the i > 1000 regime, because the beam uncertainty is smaller than for the 
first flight. The 150 GHz detectors were polarization sensitive bolometers of the 
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sort that will be used on Planck, and polarizing grids were used with the 240 and 
340 GHz bolometers: good multi-band EE and TE spectra should emerge as well. 

The South-Pole-based DASI (the Degree Angular Scale Interferometer) has 13 
dishes of diameter 0.2 m and uses high-electron-mobility transistors (HEMTs) span- 
ning 26-36GHz. An interferometer baseline directly translates into a Fourier mode 
on the sky. The dish spacing and operating frequency dictate the I range. In DASI's 
case, the range covered is 125 < I < 900. The total area covered was 288 deg 2 con- 
sisting of 32 independent maps of size 3.4°, the field-of-view. In April 2001, DASI 
unveiled a spectrum (Halvcrson et al. 2002) similar to that reported by Boomerang 
at the same time. The two results together reinforced each other and lent consid- 
erable confidence to the emerging (January 2002) d spectrum in the I < 1000 
regime. 

Technical notes 

For Maxima, 13 bandpowers, covering 36 < t < 1235, were used, with a 4% 
calibration uncertainty and a 5% beam uncertainty. For both Boomerang sets, the 
power spectra use the 'Faster' method on 3.5' pixels, the calibration error is 10% 
and the beam uncertainty is 13%. Of the 25 bands, 2-21 were used in our parameter 
analyses, covering 26 < I < 1025. The other bands are marginalized. For DASI, 
the calibration uncertainty is 4%, there is no beam uncertainty, and all 9 bands are 
used in our analyses. 

(c) June 2002: CBI, VSA and BIMA 

The CBI (Cosmic Background Imager) is based 17000 feet above sea level on the 
Atacama Plateau in Chile. It has 13 dishes, 0.9 m in diameter, operating in the same 
HEMT channels as DASI. The instrument measures 78 baselines simultaneously. 
The larger dishes and baselines imply higher resolution, to t of 3500, a huge increase 
over Boomerang, Maxima and DASI. Only the analyses of data from the year 
2000 observing campaign were reported in May 2002, covering three deep fields, of 
diameter ca. 0.75° (Mason et al. 2003), and three mosaic regions, each of size ca. 13 
deg 2 (Pearson et al. 2003). Mosaics lace together multiple interferometer fields-of- 
view (FOVs) which overlap, allowing greatly improved resolution in A£ over what 
many independent FOVs can give. 

CBI data from 2001 roughly doubled the total amount available and increased 
the area covered. The combined two-year data are for three mosaic fields covering 
about 80 deg 2 , including the two deep fields which overlap, and for one disjoint 
deep field. Preliminary results for the CBI 2000+2001 data were reported at the 
Royal Society meeting (January 2003). A full analysis is given in Rcadhcad et al. 
(2003). 

The Very Small Array (VSA) in Tencrife, also an interferometer operating at 
30 GHz, covered the I range of DASI, and confirmed the spectrum emerging from 
the Boomerang, Maxima and DASI data in that region in May 2002 (Scott et al. 
2002). It is part of the June 2002 set. The VSA observed at longer baselines to 
increase its t range and announced results in December 2002 (Grainge et al. 2003). 
It is included in the January 2003 set. 

In June 2002, the interferometric millimetre array BIMA, operating instead at 
30 GHz, announced power at I ~ 6000 had been found in random fields (Dawson et 
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al. 2002). This is too far out in t to be used for the analysis of primary anisotropies, 
but is important for secondary anisotropies (see § 2g). 

Technical notes 

The one-year CBI mosaic and deep field data use a conservative calibration 
error of 5%. For parameters, bands 1-12 covering i to 1940 are used in the odd 
binning of the mosaic power spectrum of Pearson et al. (2003). The 10 bands at 
higher £ are marginalized. When the deep field data are used in conjunction with 
the mosaic for optimal spectra, all bands are used. The deep data are not used for 
cosmic parameter estimation. The two-year data had a 3.3% calibration error as 
of January 2003, but with WMAP's observation of Jupiter this could be reduced 
to 1.3%, with a 3% recalibration downward in amplitude. The June 2002 VSA 
results have a calibration uncertainty of 3.5%, and have 10 bands up to I of 1000. 
The January 2003 extended VSA results have the same calibration, and 16 bands 
extending to I of 1700. For March 2003 they were recalibrated with the lower 
calibration uncertainty, as for CBI. 

(d) January 2003: ARCHEOPS and ACBAR; Extended Boomerang, VSA & CBI 

The January 2003 data mix replaces Boomerang, the VSA and CBI by their 
extended sky results described in § 2b, c. 

ARCHEOPS mapped 30% of the sky at a FWHM resolution of 15' in a 22 hour 
balloon flight. It used the bolometer detectors and scanning strategy similar to those 
that will be employed for the Planck satellite. In October 2002, the ARCHEOPS 
team presented a power spectrum to I = 350 that overlapped with the COBE/DMR 
power spectrum at low t (Benoit et al. 2003), derived from 12% of the sky, one of 
their six 143 GHz bolometers and one of their six 217 GHz bolometers. 

ACBAR was installed in January 2001 on the 2.1m diameter Viper telescope at 
the South Pole. It has 16 bolometers in three frequency bands centered at 150, 220 
and 280 GHz, with a best resolution of ca. 4', allowing excellent coverage to high 
I. The analysis to obtain the power spectrum of Kuo et al. (2002) used the four 
bolometers at 150 GHz available in the first observing season, the eight available in 
the second, and covered ~ 24 deg 2 in two fields consisting of correlated ~ 3 x 1.5 
deg 2 patches. The cosmological parameters determined with its power spectrum in 
conjunction with other data were given in Goldstein et al. (2002). ACBAR is being 
further upgraded and continues to operate. 

In Sept 2002, DASI announced the detection of polarization, and this is discussed 
in § 2f and as Pillar 6. 

Technical notes 

The ACBAR set has a calibration uncertainty of 10% and an uncertainty of 
3% on its <~ 5.5' beam at 150 GHz. For the analyses here, bands 2-14 are used, 
covering 300 < t < 3000, with the first bin marginalized. (The window function 
of the first bin is oscillatory and reaches reasonably high I. However, for optimal 
spectra and parameter results in Goldstein et al. (2003), the first bin has also been 
used, and yields similar results to those obtained without it, as described below.) 
ARCHEOPS has a calibration uncertainty of 7% and an uncertainty of 10% on its 
15' beam. Sixteen bands covering 15 < I < 350 are used. 
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(e) March 2003: WMAP 

In June 2001, NASA launched the five-frequency all-sky HEMT-based WMAP 
satellite. The best resolution of 12.6' at 94 GHz corresponds to a 'Gaussian' beam 
size of £ s=s 640. The first-year results unveiled in February 2003 were as spectacular 
as forecast, definitive through the second peak. The noisy error bars at higher £ 
will subside as the observing period increases, and the beam is very well known so 
that we may expect good spectra out to £ <~ 1000 with four years of data. 

COBE was fundamental to every parameter determination with pre- WMAP 
data. The WMAP verification that the COBE/DMR maps were accurate in detail 
was by itself an important step. Pillar l's Sachs- Wolfe effect dominates at low £. It 
includes both the 'ordinary' effect from fluctuations in <!> on the last scattering sur- 
face and the 'integrated' Sachs- Wolfe (ISW) effect, from the change of $ with time 
so the energy of a photon climbing out of a potential well differs from that when it 
dropped in. Other effects can influence low £, the contribution from tensor modes, 
and more exotic possibilities involving low mass scalar fields, modified topologies, 
radically broken scale invariance, etc. For the best-fit ACDM models, an upturn in 
Ct from the ISW is predicted, but instead a downturn is observed with WMAP, con- 
firming and extending the puzzles associated with the relatively low quadrupole. It 
has been our practice to drop the £ = 2 mode in parameter estimation (marginalize 
over it), and we continue it here, but recognize it could be pointing to new physics. 
The quadrupole has such large cosmic variance that including it or not does not 
change parameter determination in the minimal inflation models much. 

The low-£ DMR data has been used to constrain the size of the Universe, basi- 
cally from the scale of its hot and cold spots. If our manifold was much smaller than 
the apparent distance to the last scattering surface, there would not be large-scale 
spots in DMR maps, and smaller scale ones would be images of each other, with 
placement and details dependent upon specifics of the manifold. The complication 
here is that there are many possible manifolds, and each has an orientation, so it 
is possible to find universes that are just big enough to be more probable than 
the conventional inflation models with far fewer degrees of freedom in which the 
large-scale hot and cold spots are not geometrically correlated as in the topology 
case. See Bond, Pogosyan & Souradeep (2000) for a discussion. WMAP data should 
allow more powerful tests and improve the current constraints on size somewhat. 

Technical notes 

The WMAP power spectra for our analyses are those obtained in Hinshaw et 
al. (2003) for TT to £ = 900 and in Kogut et al. (2003) for TE to £ = 512, using 
the construction in Verde et al. (2003) for the correlated errors. 

For our Cg database approach, we include the 0.5% calibration uncertainty, use 
the offset lognormal approximation to the likelihood surface, and, as for DMR, we 
marginalize over the WMAP quadrupole. We have compressed the data to speed 
up our analyses. We have used the full 899 individual multipoles of Hinshaw et al. 
(2003) and Verde et al. (2003), including correlations out to a (tiny) sideband cutoff. 
Doing many large matrix multiplies can be slow, as relaxation to the amplitudes 
is done for each of the 8.5 million model elements in the database. Although these 
operations could be sped up considerably by approximate log-likelihood estimates 
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to reject extremely improbable models, and there are many of these, we have so far 
used brute force through all models. 

We have also used compressions of the WMAP data onto fewer bands for pa- 
rameter estimations, optimal spectra and calibration estimation. For example, a 
98-band compression of the data preserves all signal and noise information and 
side-band correlations in the bands, has A£ spacings ranging from 1 at small £, to 
7 to 5 to 7 through £=373, 9 through 565, with a gradual increase beyond 600, 
and only 4 broad bands beyond 700 (with relative errors above 50%). A best fit 
shape was used in the compression, but other choices, including a flat shape, give 
very similar spectra where WMAP dominates the data. In a 49-band compres- 
sion, the band spacings are essentially doubled, and in a 19-band compression, the 
Boomerang A£ = 50 spacing was used. Parameter results using the Cg database 
are quite consistent among the 899-band, 98-band and 49-band sets, with slightly 
sharper errors with the 899-band set. We do find small (sub-ltr) deviations in pa- 
rameters that are bundled into some of the exquisitely determined eigenmodes for 
the 899 set compared with the 98 set. 

For our MCMC calculations, the WMAP likelihood routines of Verde et al. 
(2003) are adopted, since these are used in the Cosmomc package (Lewis & Bridle 
2002). This uses a hybrid of an offset lognormal and a Gaussian distribution to 
compute parameters. The quadrupole is included. The small calibration error of 
0.5% is not. For Cg construction, CMBfast was used for the database and CAMB is 
used in Cosmomc. Other differences between the database and MCMC treatments 
of WMAP are bigger, e.g. mimicking the effects of the TE data with a tq prior. 
The parameter results for the 98-band case using MCMC are in good agreement 
with those using the database. 



(/) Beyond March 2003: Planck and targeting polarization 

Many other CMB experiments on the ground and in balloons will happen before 
the European Space Agency's (ESA's) Planck satellite is launched in 2007. Planck 
combines bolometers, many of which will be polarization sensitive, and HEMTs. Its 
best resolution of ~ 5' and detector sensitivity should allow the damping tail and 
power spectrum modifications due to secondary effects, including weak lensing, to be 
very well determined. The all-sky component-separated intensity and polarization 
maps will be superb for interstellar matter and extragalactic source research as well 
as cosmological research. Such accuracy is needed to open up the cosmic parameter 
space and search for anomalies that may signal new physics beyond the minimal 
inflation concordance model we are drawn to now. 

Polarization is described by a 2x2 tensor on the sky, with the components 
related to the four Stokes parameters, T,Q,U,V, with T the total intensity (i.e. 
temperature) . The polarization dependence of Compton scattering induces a well- 
defined polarization signal emerging from photon decoupling, arising from the 
quadrupole nature of the viscosity-induced anisotropic stress tensor. For primary 
CMB fluctuations, circular polarization is not there and V vanishes. When Q and 
U maps are Fourier transformed and are rotated into a basis related to the angular 
wavevector, these give £7- type (grad) and £?-type (curl) 'maps'. Scalar perturba- 
tions in linear theory do not generate £>-type modes, so their absence is a check. 
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Tensor perturbations generate both, and detection of B-modes at low I would be 
a direct signature of a gravitational wave background (pillar 7). 

Given the total C^ TT ^ of figure 1, we can forecast the polarization power c[ EE ^ 

and cross-correlation power Cf E) : the maximum signal is expected at I <~ 900, 
with amplitude <~ 5^iK over t <~ 400 — 1600. 

The great race to first detect CMB E-m.oAc polarization was won by DASI 
(Leitch et al.2002, Kovac et al. 2002), with 271 days of polarization data on rwo 
deep fields (3.4° FOV) showing a 5a detection with a value 0.8 ± 0.3 of the fore- 
casted amplitude from T for inflation-based models. The cross-correlation of the 
polarization with the total anisotropy had an amplitude 0.9 ± 0.4 of the forecast. 
These detections used a broadband shape covering the I range <~ 250 — 750 de- 
rived from the theoretical forecasts. The powerful cosmological implications of the 
remarkable WMAP TE bandpowers have already been discussed. 

Forecasts indicate solid EE power spectrum determinations are likely soon from 
the ongoing CBI polarization observations and Boomerang's January 2003 flight 
with polarization-sensitive bolometers. Both are optimally sensitive to the t <~ 900 
region where the EE power is expected to peak. MAXIMA will fly again as the 
polarization-targeting MAXIPOL. Other EE experiments, operating or planned, in- 
clude AMiBA, CAPMAP/PIQUE, COMPASS, CUPMAP, POLAR, Polarbear, Po- 
latron, QUEST and Sport/BaRSport, among others. We are also awaiting WMAP's 
EE results. 

Although the strength of the S-mode induced by gravity waves is model specific, 
the amplitude is expected to be quite small even at low I. Nonetheless there are 

(BB) 

experiments such as BICEP being planned to go after C\ in these low I ranges. 
Planck could also make such a detection. A nice figure summarizing EE and BB 
bandpower forecasts for various experiments is given in Hivon & Kamionkowski 
(2002). The promise is sufficiently exciting that a CMBPol satellite is being con- 
templated by NASA as the next step in space for the CMB after Planck. 

(g) Beyond March 2003: targeting secondary anisotropics 

Spectral distortions from the CMB black body must exist as a result of nonlin- 
ear processes and will have associated anisotropics. The spectrally well-defined SZ 
distortion associated with Compton upscattering of CMB photons from hot gas has 
not been observed in the spectrum. The COBE/FIRAS 95% upper limit of 6 x 10~ 5 
of the energy in the CMB is compatible with values up to around 10 -5 expected 
from the cosmic web of clusters, groups and filaments in the inflation-based class 
of models considered here, and places strong constraints on the allowed amount of 
earlier energy injection, e.g. ruling out mostly hydrodynamic models of LSS. 

The SZ effect has been well observed at high resolution with very high signal-to- 
noise along lincs-of-sight through a large number of clusters now, by single dishes, 
the OVRO and BIMA millimetre arrays, and the Ryle interferometer. This tells us, 
among other things, that the CMB comes from further away than redshift z <~ 1 
- if we had any residual doubt. The SZ effect in random fields may be responsible 
for the power at I > 2000 seen in the CBI deep data (figure 3) (Mason et al. 
2003; Bond et al. 2003b), in the BIMA data at I ~ 6000 (Dawson et al. 2002), and 
possibly in the ACBAR data (Goldstein et al. 2003). Multi-frequency observations 
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to differentiate the signal from the CMB primary and radio source contributions 
will be needed to show this. 

A number of planned HEMT-based interferometers are being built with this am- 
bient effect as a target: CARMA (OVRO+BIMA together), the SZA (from Chicago, 
and to be incorporated in CARMA) , AMI (based in Britain, including the Ryle tele- 
scope) and AMiBA (from Taiwan). Bolometer-based experiments will also be used 
to probe the SZ effect, including: the CSO (Caltech Sub-Millimetre Observatory, a 
10 m diameter dish) with BOLOCAM on Mauna Kea; the LMT (Large Millimetre 
Telescope, with a 50 m diameter dish) in Mexico; APEX, a 12 m diameter German 
single dish based in Atacama; Kobyama, a 10 m diameter Japanese single dish; 
and the 100 m Green Bank telescope. Large bolometer arrays with thousands of 
elements and resolution below 2' are also under development: the South Pole Tele- 
scope (SPT, Chicago) and the Atacama Cosmology Telescope (ACT, Princeton). 

The kinetic SZ effect due to the motion of clusters in the low-to-moderate- 
rcdshift cosmic web, or just clumpy ionized gas in the high redshift cosmic web, has 
the same spectral signature as primary anisotropies, and so is harder to disentangle. 

The (non-Galactic) distortions from the black body that have been detected in 
the COBE FIRAS and DIRBE data are associated with starbursting galaxies due 
to stellar and accretion-disk radiation downshifted into the infrared by dust then 
rcdshifted into the sub-millimetre. This background has energy about twice the 
total of that in optical light, about a tenth of a percent of that in the CMB. About 
50% of this sub-millimetre background has been identified with sources found with 
the SCUBA bolometer array on the JCMT. Anisotropies from dust emission from 
these high redshift galaxies are being targeted by the JCMT (with a SCUBA-2 
very large array in development), the OVRO millimetre interferometer, the CSO, 
the SMA (Sub-Millimetre Array) on Mauna Kea, the LMT, the ambitious US/ESO 
ALMA millimetre array in Chile, the LDB BLAST, and ESA's Herschel satellite. 

3. CMB analysis and phenomenology 

(a) CMB pipelines: from timestreams or visibilities to bandpowers 

Analysing CMB experiments involves a pipeline that takes the raw detector 
timestreams or visibilities from correlators for interferometers, flags and cleans 
them, and usually generates maps, from which bandpowers and higher order statis- 
tics are derived, ideally after separating component signals by using their differing 
frequency and spatial dependences. The step from bandpowers to cosmic param- 
eters described in § lc may be the goal, but it is not where most of the time is 
spent. Indeed parameter determination is used as a diagnostic along with every- 
thing else as the CMB teams struggle to understand in detail their experimental 
results. Every new round of data generates a fresh look at pipelines, and often 
new faster algorithms for proceeding. Recent pipelines are described in e.g. Bond 
& Crittenden (2001), Netterfield et al. (2002), Hivon et al. (2002) and Ruhl et al. 
(2003) for single-dish bolometers, in e.g. Myers et al. (2003) for interferometry, and 
in Bennett et al. (2003), Hinshaw et al. (2003) and Verde et al. (2003) for WMAP 

For single-dish experiments, the timestreams are turned into spatial maps for 
each frequency: an average temperature in each pixel and a pixel-pixel noise corre- 
lation matrix from which the bandpowers, noise in the bandpowers and band-band 
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error matrices are derived. The first step is to extract the sky signal from the noise, 
using the only information we have, the pointing matrix mapping a bit in time onto 
a pixel position on the sky. In the analysis of Boomerang, and subsequent work in- 
cluding for WMAP, powerful use of Monte Carlo simulations was made to evaluate 
the power spectrum and other statistical indicators in maps with many more pixels 
than was possible with the conventional matrix methods described in Bond, Jaffe 
& Knox (1998) and Borrill (1999). 

For interferometer experiments, the basic data are visibilities as a function of 
baseline and frequency, with contributions from random detector noise as well as 
from the sky signals. A baseline is a direct probe of a given angular wavenumbcr 
vector on the sky, and hence suggests that we should make 'generalized pixel' maps 
in 'momentum space' (i.e. Fourier transform space) rather than in position space, 
as for Boomerang. For CBI, the > O(10 5 ) visibility measurements of each field were 
'optimally' compressed into a < 0(1O 4 ) coarse-grained lattice in momentum space 
from which the power spectrum was calculated using matrix methods (Myers et al. 
2003). 

The important step of separating multi-frequency timestream data into the 
physical components on the sky is fundamental, still under active development, 
and will remain so for a long time, as our precision increases. The sources are 
the primary CMB, the thermal and kinetic SZ effects, the dust, synchrotron and 
bremsstrahlung Galactic signals, the extragalactic radio and sub- millimetre sources, 
possibly spinning dust, and of course the sources we have not thought of yet. An 
example is the treatment of point sources: at the high resolution of CBI and its 30 
GHz frequency, the contribution from extragalactic radio sources is significant, so 
known point sources were projected out of the visibilities by using a large number 
of constraint-template matrices which marginalize over all affected modes, using 
positions from the (1.4 GHz) National Radio Observatory Very Large Array Sky 
Survey (NVSS) catalog (Mason et al. 2003). Even though the WMAP resolution 
is substantially lower than that of CBI, WMAP cut out some 700 sources in its 
analysis. Boomerang cut out three quasars. 

For extended sources, the most expedient method is to just cut out problematic 
areas, e.g. with Galactic latitude cuts to remove emission from the plane, as was 
done for DMR, WMAP and Boomerang. Extensive use is also made of templates 
that come from Galactic observations, e.g. the IRAS/DIRBE maps of sub- millimetre 
emission, HI maps, etc., and from the experiments themselves, with the highest and 
lowest frequency maps often being completely dominated by foregrounds (as for the 
dust-dominated 400 GHz Boomerang map). Having errors on primary signal maps 
that reflect the residual contamination after separation is a priority to ensure that 
the precision is unbiased as well as high. 

The extension of these methods to polarization is under active development. 
For CBI, the polarization pipeline using generalizations of the Myers et al. (2003) 
techniques has been implemented. For Boomerang, the basic tools being developed 
include a generalization of the signal-noise separator used in Boomerang-98 (Prunct 
et al. 2001; Dore et al. 2002), and a variety of power spectrum estimators such as 
Spice (correlation function techniques (Chon et al. 2003)), Faster (approximate t- 
space techniques (Netterfield et al. 2002; Ruhl et al. 2003)) and Madcap (full matrix 
methods (Borrill 1999)). 
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(6) Calibrating by power spectra 

The maximum-likelihood values of the calibration and beam parameters and the 
errors on them that are determined in the calculation of the optimal power spectra 
of § la are themselves very useful: for the power spectra to be consistent with one 
derived from an underlying Gaussian model for the anisotropies, the experimentally 
quoted uncertainties are often reduced, as shown in table 1. Even without the 
excellent overall amplitude accuracy to 0.5% of WMAP, self-calibration among the 
pre-WMAP experiments was giving some of the same factors, which is why the 
grand unified spectra look so good relative to WMAP. Overall figure 2 shows an 
amazing concordance of data. Direct calibration using source observations or maps 
on equivalent regions is preferable. This has been done for Boomerang with WMAP 
maps by Hivon (w 0.95 ± 0.02) and for CBI with WMAP observations of Jupiter, 
the CBI calibrating source. It is remarkable how close these recalibrations are to 
those given here using global power spectrum analyses. 

The CBI one-year data started with a ±0.05 calibration error. Running the 
optimal spectrum with WMAP, CBI is predicted to have its calibration lowered by 
0.970±.037. For the two-year CBI data, we begin with the January 2003 uncertainty 
of ±0.033 and get 0.970 ± .022. The Jupiter calibration with WMAP resulted in a 
3% downward calibration with a 1.3% calibration error (Readhead et al. 2003). The 
same calibration adjustment and uncertainty apply to the VSA, which also used 
Jupiter. 

Technical notes 

The WMAP power spectrum used for this calibration and beam analysis was 
the 49-band compression of the i < 900 spectrum and the optimal compression 
was onto the 18 bands of June 2002 in figure 2, but the results are very robust: 
other compressions for either WMAP or the optimal spectra, varying the Ci shape of 
figure 2 used to get the bandpowers, and adding more experiments in the calibrating 
mix make little difference. 

(c) The phenomenology of peaks and dips 

The emerging structure in the evolving power spectra of figure 2 has stimu- 
lated a number of 'model-independent' approaches to determining the statistical 
significance of peaks and dips. Here we apply a simple procedure we have used on 
Boomerang and CBI. We model the local band-power profile as the quadratic form 

Cb = C p k + K pk ((£ — l p k) 2 )b/2 , 

where b = 1,..,N runs over the number of bands N we include as we slide over the 
data, and ()b denotes band-average. The position £ p k, amplitude C p k, and curvature 
K p k, are treated as 'internal' parameters, just as calibration and beam uncertainties 
are: they relax to their maximum-likelihood values and errors are estimated from 
the curvature matrix of the likelihood function about this maximum. A numeri- 
cal indication of the significance of a detection is the number of sigma the peak 
curvature differs from zero, \K p k. m \/o' Kpk . To detect a peak, we require that this 
significance exceeds unity, and also that £ p k, m lies within the range of multipolcs 
covered by the bands. The results of applying this algorithm to the evolving grand 
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unified spectra and to WMAP are shown in table 2. No attempt was made to 
optimize these detections by exploring band positioning and spacings. 

The first peak was hinted at in the April 1999 data, was seen in the January 
2000 data when Toco and Boomerang-NA were added, and localization improved 
by January 2002 with Boomerang, Maxima, and DASI, and then with ARCHEOPS 
by January 2003, and now to very high accuracy by WMAP. Boomerang+DASI 
(2001) detected the second and third peaks, first and second dips, and CBI (2002) 
detected the second, third, and fourth peaks and the third and fourth dips, albeit 
with some at low significance. 

Technical notes 

The slider bandwidth TV we choose depends upon the band-spacing. For the 
unified spectra of the evolving data, we compressed onto the A£=50 bins used for 
Boomerang up to 900, then onto a 'CBI odd' binning used for the two-year data, 
although the CBI bands used were for the 'CBI even' binning: the odd binning 
is more similar to the ACBAR spacing. For the June 2002 data, CBI one-year 
A£ = 140 odd bins were used. Sliders of width N — 3,4,5 were tried, N = 3 
being more appropriate for high I values and N — 4 being more appropriate for 
lower values, to keep the the Grange probed limited to a few hundred. WMAP data 
were first compressed onto this Boomerang binning before the optimal spectra were 
constructed for the March 2003 data, but results are relatively insensitive to this 
and to band locations and widths. For ARCHEOPS, the A£ are smaller, so larger 
N gives the most robust detections. 

Table 3 is another approach to peak/dip detection: given a class of theoretical 
models with a sequence of peaks and dips, the statistical distribution of positions 
and amplitudes can be predicted by ensemble-averaging over the full probability, 
the multidimensional likelihood. 

(d) Characteristic scales: sound crossing, peaks and dips, and damping 

A strong first peak followed by a sequence of smaller peaks diminished by damp- 
ing in the Ci spectrum was a long-standing prediction of adiabatic models (pillar 2). 
The critical scale determining the spatial positions of the acoustic peaks in the spec- 
tra of figure 1 is the (comoving) sound crossing distance at recombination, r s . The 
corresponding multipole scale is £ s = TZdec/fs, where IZdec is the angular-diameter 
distance that maps angles observed at our location to comoving spatial scales at 
recombination. 

In terms of the comoving distance Xdec to photon decoupling (recombination, 
at redshift Zd ec = oJ^ ec ~ -0' an< ^ the curvature scale dk, TZdec is given by 



The three cases are for negative, zero and positive mean curvature. Thus the map- 
ping depends upon oj k , uq and wq as well as on ui m . 
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The sound crossing distance at recombination is 



6000 



V3 



Mpc 




V 172 ^ + w 6 a/(4o; 7 /3))- 1 / 2 d^a~, (3.2) 



where lo 1 = 2.46 x 10~ 5 is the photon density and oj er — 1.68w 7 for three species 
of massless neutrinos. 

The estimates of r s and n£ s have been quite stable over time (table 3). Values 
are determined by averaging over the Ci model space probabilities. Since the r s are 
comoving, the physical sound horizon at decoupling is ca. 140 kpc. 

The angular-diameter-distance relation maps spatial structure at photon decou- 
pling perpendicular to the line-of-sight with transverse wavenumber k± to angular 
structure, through I = lZd ec kj_. Converting peaks in fc-space into peaks in ^-space 
is complicated by three-dimensional to 2D projection effects over the finite width of 
decoupling and also by the influence of sources other than sound oscillations such 
as Doppler terms. 

The peak locations £ p kj in table 3 are obtained by forming cxp < ln£ p kj >, 
where the average and variance of \n£ p kj are determined by integrating over the 
probability-weighted Ci database. These peak locations accord reasonably well with 
£pk,j ~ jfj 7T £ s , where the numerically estimated constant fj w 0.75 for the first 
peak, approaching unity for higher ones. (There are small n s — 1 corrections to 
hi/j.) The interleaving dips are also shown. Dip locations are well determined by 
replacing j by j + 1/2, with slightly different fj factors. 

For fixed u>b and u m , constant £ s lines in the flk~^A and wq-£Iq planes look 
rather similar to contour lines determined from the data (see figures 5 and 7 in 
Bond et al. 2003a). This degeneracy (Efstathiou & Bond 1999) among these pa- 
rameters would be exact except for the integrated Sachs- Wolfe effect. However, lines 
of constant Hq are different in that space and break the degeneracy, only weakly 
for the weak prior but more so when the HST-h or the SN1 prior on deceleration is 
imposed. The degeneracy in the f2 fc -^A plane is also broken when LSS information 
is added, though less so in the wq-^Iq plane. 

Also evident in the spectra in figure 1 and the diminishment of the peak heights 
in table 3 is the damping tail, an overall decline due to the shear viscosity and 
the finite width of the region over which hydrogen recombination occurs. Scales 
for both can be estimated analytically and are similar. Values for the comoving 
scale Rd and an associated angular damping scale Id = TZdec/RD are given in 
table 3. The physical scale Rd/(1 + Zd ec ) at recombination is therefore about 10 
kpc. (The decline in Cg is parameterized by a exp[— [£/£D) mD ] multiplier (Hu and 



White 1996). We find m D = 1.27±g;ggl for January 2003 and 1.27±g;gg| for March 



The dependence in r s would lead to a degeneracy with other parameters 
in terms of peak/dip positions. However, relative peak/dip heights are extremely 
significant for parameter estimation as well, and this breaks the degeneracy. For 
example, increasing ujb beyond the nucleosynthesis (and CMB) estimate leads to a 
diminished height for the second peak that is not in accord with the data. Page et 
al. (2003) used the peak/dip parameters to estimate directly the cosmic parameters. 



2003.) 
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4. Cosmic parameter estimations 

(a) Evolution of marginalized cosmic parameters 

Table 4 and figure 4 show the evolution of cosmic parameter estimations from 
projected ID likelihood curves, using the Ci database used in Sievers et al. (2003), 
Ruhl et al. (2003) and Goldstein et al. (2003), an extension of that used in Lange et 
al. (2001), JafFe et al. (2001) and Netterfield et al. (2002) to reflect the growing pre- 
cision. The ranges quoted are Bayesian 50% values and the errors are la, obtained 
after projecting (marginalizing) over all other parameters. For derived variables 
such as h amd age, ensemble-averaged means and variances are used. The WMAP 
addition punches out the detection in a remarkable way. 

The MCMC parameter estimates for weak and weak+flat priors in table 5 agree 
well with the entries in table 4: about as good as one might expect given the dif- 
ferences. In particular, WMAP with MCMC includes the tq constraint by explicit 
calculation of TE likelihoods, while the database either lets the TT data alone de- 
cide or uses the broader-than-Gaussian tc prior to reflect the detection. The March 
2003 results in the table include this prior, but most parameters are insensitive to 
it: even cr 8 does not migrate that much. 

The Spcrgel et al. (2003) MCMC four-chain results for flat power law ACDM 
models for WMAP-only are quite similar to what we obtain. Their WMAP-ext 
consists of WMAP plus cut verions of the ACB AR and CBI one-year power spectra. 
Our March 2003 set contains more data and allows more overlap. Their analogue of 
LSS is to use the 2dFRS power spectrum. This assumes a linear scale-independent 
bias. It is a more stringent prior on shape than our T constraint. They also add 
(small-scale) Lya forest information to extend the fc-space coverage. The forest 
data could be sensitive to gastrophysical complications, and we would hesitate to 
construct any but a very weak prior for it at this time. 

We also find good agreement with Contaldi, Hoekstra and Lewis (2003), who 
applied Cosmomc to the WMAP+ACBAR+CBI data and the Red Cluster Survey 
(RCS) weak-lensing data. Our LSS prior was designed to encompass as constraints 
from weak lensing displayed in figure 5, and RCS is one of those entries. 

Table 4 shows that with just the weak prior, there are strong detections for Sltot, 
Wb, ^cdm and n s . The HST-h prior helps to determine better than the weak prior, 
because it breaks more strongly the flk^A near-degeneracy. SN1 breaks it even 
more strongly, and the HST-h, weak+SNl and weak+LSS results are all compatible, 
showing there is no very strong tension among these priors. 

The precision of WMAP invites exploration of larger parameter spaces. A first 
issue is whether n s varies. When Spergel et al. (2003) added the high-£ ACBAR 
and CBI data to make their WMAP-ext dataset, they found lower n s , and further 
adding the 2dFRS data and the small-scale data from the Lyman alpha forest 
exacerbated the issue. We find the same result with the database with MCMC: n s 
changes by about 0.03 from WMAP-only to the March 2003 set for the weak+flat 
prior, and by a further 0.015 with LSS added. 

Modelling n s variation with a logarithmic correction, Spergel et al. (2003) get 
n s varying from 1.2 to 0.93. (See also Bridle et al. 2003, who show that dropping 
low-£ modes reduces the significance of running index detection.) The lower panel of 
table 5 shows how the parameters change relative to the fixed n s model when this ex- 
tra parameter is introduced, for the weak+flat prior and for the various datasets. We 
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restricted dn s /dh\k to lie within the range —0.2-0.2. The normalization point for 
n s is k- 1 = 20 tr 1 Mpc, and the horizon scale is \dec « 0.88(6000/fi^ 2 ) h" 1 Mpc 
for a ACDM model, so In k n Xdec ~ 6, hence substantial slope changes are possible 
for the dn s /d\nk values obtained. Although table 5 shows there is a preference for 
a logarithmically varying slope from January 2002, due to downturns in the CMB 
data at low and high I, it is < 2a in likelihood fall-off even in the March 2003 data. 
The likelihood maxima are also smaller than the means given in the table. We 
also caution that dn s jd In k(k n ) is correlated with the Tc-o-$-n s combination, which 
results in upward extension of the ID distributions for tc and as- Prior choices 
on tc can then move the limits around: the weak+flat results in table 5 place no 
restriction beyond the weak 0.7 upper boundary on tc to conform MCMC to our 
database limit. 

The WMAP team estimate the gravitational wave (tensor) contribution to be 
< 0.72 of the scalar component in amplitude. 

We have considered variations in the dark energy equation of state. For the 
database calculations, we used a flat prior and allowed wq to vary from —1 to 
—0.01, as in Bond et al. (2000). We applied it to the June 2002 data (and obtained 
wq < -0.7 at the 95% CL for the prior combination weak+flat+LSS+SNl) (Bond 
et al. 2003a) and to what was almost the same as the January 2003 data (for which 
wq < —0.8 was obtained) (Pogosyan et al. 2003). MCMC calculations yield about 
the same limits: less than —0.70 at 2a for the HST-h+flat+LSS and SN1 prior. 
We get only slightly better limits for the March 2003 data (less than —0.71). This 
constraint is SNl-driven (for weak+SNl wq < —0.68) rather than HST-h-driven 
(wq < —0.48). Although bands of allowed h do break the angular-diameter-distance 
degeneracy between Qq and wq, allowed bands for the deceleration parameter q 
break it more effectively. 

Unlike the WMAP team, our weak+flat+LSS prior does not give a good con- 
straint. Since the shape parameter L <~ u m /h, and u) m is now accurately determined 
by the CMB data, our shape constraint approaches a pure h prior. For w m <~ 0.14, 
it is similar to the weak limits we are already imposing on h. The 2dFRS data give 
T = 0.21 ± 0.03 and the SDSS data give 0.19 ± 0.04, so the combination would 
nomimally have an 11% relative Gaussian error, the same as for the HST-h prior, 
suggesting the constraint would be similar if this strong-LSS were to be imposed. 

Why the limits are typically < —0.7 or —0.8 with SN1 is easy to understand: 
the 2a SN1 error contour in the £Iq-wq plane roughly follows a q ~ — 1/4 line, 
where qo = (1 + 3wq£Iq)/2 is the current deceleration parameter. With fig > 0.65, 
this gives the -0.7 to -0.8 range for the limit. The way the CMB comes into this is 
to restrict the allowed values of Qq , and it used to be that the data was such that 
the LSS prior was needed to do this the first time we went through this exercise. 
The addition of WMAP refines Qq but this only fine tunes the limit by a small 
amount. 

The other component to our LSS prior, the amplitude constraint as, has had 
a significant impact in constraining models. The June 2002 data for the weak + 
flat prior give a s = 0.88 ± 0.11 and when LSS is added, 0.86 ± 0.08 (Bond et 
al. 2003b). For the March 2003 data, the MCMC calculation with the wcak+flat 
prior gives 0.85 ± 0.06, the database calculation with the weak+flat+Tc prior gives 
0.83 ± 0.06 using the 899-point WMAP data and 0.81 ± 0.06 using the compressed 
98-band WMAP data. Adding the LSS prior gives 0.86±0.05, 0.85±0.04 and 0.86± 
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0.04, respectively. (These a s differences, plus ~ 0.01 excursions in n s , represent the 
biggest deviations between compressed and uncompressed WMAP estimations; and 
the up and down la and 2a limits are close.) 

The SZ effect breaks the a%-T C degeneracy. The SZ power spectrum is found 
to scale as a\ about a ACDM model with a s = 0.9 (e.g. Bond et al. 2003b). This 
high non-linearity means that it could be an excellent way to estimate as- For an 
SZ explanation to work for the (recalibrated) CBI+BIMA+ACBAR data seems 
to require as ~ 0.94lg;°g with la errors which include the non-Gaussian nature 
of the SZ effect in the small patches that the CBI deep and BIMA measurements 
probe (Goldstein et al. 2003). These SZ distributions overlap with those found for 
the CMB and CMB+LSS data, but the jury is still out on whether this is the 
explanation. 

The calculation of the SZ angular power spectrum requires detailed hydrody- 
namical simulations that properly take into account cluster structure, pressure pro- 
files, heating/cooling, etc. There are still uncertainties as to the role of feedback, 
and much work is needed on the theoretical side as well as the observational side 
to turn this technique into a high-precision tool. 

(b) Parameter eigenmodes and degeneracy breaking 

To help break parameter near-degeneracies via CMB alone, one can decrease A£ 
of the bands with fixed error, as ARCHEOPS and WMAP did. We can also extend 
the £ range, which CBI, ACBAR and the VSA have done. Planck will do both. 

Polarization helps, by breaking the r^-amplitude near-degeneracy and also by 
breaking the n s (k) degeneracy that arises if we let the index have full functional 
freedom. The EE power near the peak will be well probed by many of the polar- 
ization experiments mentioned in § 2f. 

Degeneracy breaking also arises when non-CMB data are added, as the HST-h, 
SNla, and LSS examples show. Higher order statistics can also help, e.g. skew- 
ness/kurtosis. 

Forecasting how proposed experiments would improve on cosmic parameter er- 
rors became, and remains, quite a cottage industry. As mentioned in § le, a figure of 
merit is how many parameter eigenmodes will be determined to a specific precision 
level, i.e. the number of modes with a a below some number, taken here to be 0.1 
and 0.01, respectively. 

How well did we do? The forecast for Boomerang+DMR for 1.8% of the sky 
with four 150 GHz bolometers gave four out of nine linear combinations should 
be determined to ±0.1 accuracy. This is what was obtained in the full analysis 
of Netterfield et al. (2002). (For the forecasts, oJhdm and n t were included in the 
parameter mix, wq was frozen at —1, making nine.) The forecast for WMAP with 
two years of data was 6/9 to ±0.1 and 3/9 to ±0.01. For one year of WMAP-only 
data and the weak prior, 2/7 are determined to ±0.01, 5/7 to ±0.1; adding the tc 
prior increase this to 6/7. (For Planck, 5/9 to ±0.01 accuracy are predicted.) 

The March 2003 eigenmodes for the weak prior are: to ±0.004, a primarily f^- 
17a combination, related to ir£ s ; to ±0.008, a \nCio-Tc-n s combination; to ±0.02, 
uib predominantly, with some extra; to ±0.04 a In C\$-Tc-n s -u) c dm combination, and 
another to ±0.08. The sixth, ±0.13, is in a direction that would break the amplitude 
tc degeneracy. When the tc prior is included, this error diminishes. The least well 
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determined mode is associated with the fi^-f^A angular-diameter-distance degener- 
acy. Goldstein et al. (2003) give the eigenmodes for what is basically the January 
2003 mix except the June 2002 CBI and VSA data are used. The January 2003 
data used here have best determined modes, to ±0.011, ±0.015, ±0.05, that are 
mainly similar to the March 2003 ones, but with differing subdominant structure 
in the linear combinations. 

Whether before or after WMAP, the simplest inflationary paradigm with min- 
imal parameters fits the data well. This docs not mean inflation is proved, but 
competitor theories would have to look awfully like inflation for them to work. As 
CMB precision increases at high £, and more polarization data arrive, and as the 
LSS and SN1 data improve, more details with more parameters will be explored. 
Already WMAP points once again to the low-£ anomaly and with the other data 
to possible variation in the primordial slope with wavenumber. 

We thank our collaborators on the Boomerang, CBI and ACBAR teams for sharing the 
exhilaration of spectral revelation over the past few years and Antony Lewis for Cosmomc 
and other discussions. Lastly we would have you picture a beaming Dave Wilkinson with 
Sunyaev and Bond in front of a Saturn V rocket of the sort that took the US to the 
Moon, taken shortly before the Delta-launch of WMAP in June 2001, and the celebration 
thereafter — still ongoing, but greatly missing the W in WMAP. 
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Figure 1. Optimal Ce spectra for the pre- WMAP January 2003 data and the post-WMAP 
March 2003 data show good agreement. These spectra are maximum-likelihood deter- 
minations of the power in 26 (top-hat) bands, with calibration and beam uncertain- 
ties of the various experiments fully taken into account. Two Ce ACDM models from 
the database are shown. The dotted (black) curve best fits the March 2003 data with 
the weak+flat+LSS-|- tc prior applied. It has parameters {fitot, ^a, Q.bh 2 , Q, c dmh 2 , n s , 
T C ,t ,h,a$} = {1.0,0.7,0.0225,0.12,0.975,0.15,13.7,0.69,0.89}. The solid green curve 
that looks quite similar best fits the June 2002 data for the weak+flat+LSS prior (Sievers 
et al. 2003), with parameters {1.0, 0.5, 0.020, 0.14, 0.925, 0, 14.4, 0.57, 0.82}. It was used as 
the inter-band shape for this optimal bandpower determination, but the results are insen- 
sitive to this. The bandpowers are optimally placed in i. Their finite horizontal extension 
is not shown, and the vertical diagonal bandpower errors also do not show the whole story 
since there are band-to-band correlations, (e.g. the visual up-down-up at the first peak 
for January 2003 is indicative of the strong correlations and can disappear with different 
banding, e.g. one better tuned to Boomerang's binning.) Despite these caveats, the best 
fit Ce would fit better with a slight downward tilt beyond £ > 500, which a scale-dependent 
n s (k) could do (see § 4). 
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Figure 2. The evolving optimal Ce values are compared with WMAP-only power spectra 
(crosses) compressed onto the same bands. The degree of compression does not do visual 
justice to WMAP because the errors are so small until beyond the second peak. The bands 
were chosen to be natural for the data at the time, but band-to-band correlations do exist. 
The spectra include the following data: January 2000 has DMR + Toco + Boomerang-NA 
+ the April 1999 mix; January 2002 has January 2000 and the Boomerang data of Net- 
terfield et al. (2001) + Maxima + DASI; June 2002 has January 2002 plus one-year 
CBI mosaic and deep field data and the (non-extended) VSA data; January 2003 has 
ARCHEOPS + ACBAR, uses the Ruhl et al. (2003) Boomerang spectrum covering 2.9% 
of the sky, the extended- VSA data and the two-year combined CBI mosaic + deep field 
data. The £ > 2000 excess found with the one-year deep CBI data is denoted by the light 
blue hatched region (95% confidence limit) in the right hand panels. The two best fit 
ACDM models of figure 1 are repeated in each of the panels. When HST-h or SN1 priors 
are included in the June 2002 data, the best fit model has the same parameters as those 
of the March 2003 curve, except for a slight shift in tilt, to n s = 1.0, a corresponding rise 
in tc, to 0.20, leading to <T8=0.91. 
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Figure 3. The evolution of 2a likelihood-contour regions, with the weak+LSS prior prob- 
ability applied in the top-left-hand panel, and the flat fit t=l applied additionally in the 
rest. The outer (red) contour is for COBE-DMR only, then the sequence of increasing con- 
centration is for the January 2000 (green), January 2002 (magenta), June 2002 (cyan) and 
January 2003 (blue) data. The dense (black) region is for March 2003, including WMAP, 
but excluding ARCHEOPS and DMR because of their substantial overlap with WMAP. A 
tc prior (see § Id) taking into account the WMAP 'model independent' TE analysis, has 
been included, but it does not make much difference to these results. If only the weak prior 
is imposed, the January 2000 data still show rough Qk localization, but do not constrain 
the other parameters significantly, whereas the January 2002 and January 2003 contour 
regions without LSS are only slightly bigger than those shown here. For March 2003, the 
black regions remain small, but the Qa~ Qk degeneracy becomes more evident, with closed 
universes and smaller Qa allowed. 
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Table 1. Calibration parameters that the temperature anisotropies should be multiplied by, 
as determined in the computation of optimal spectra. (Before WMAP, calibrations by this 
method were made relative to each other, with the interferometry data (INT) dominating 
because of smaller uncertainties: 4% for DASI; 5% for CBI in the June 2002 data, 3.3% 
in January 2003, and 1.3% in March 2003 (Readhead et al. 2003); 3.5% for VSA for June 
2002 and January 2003, and 1.3% in March 2003. WMAP gives 0.5%. INTc denotes the 
recalibrated CBI and VSA, along with DASI. The Boomerang calibrations shown are for 
the 2.9% sky-cut Ruhl et al. (2003) 'Faster' power spectrum, in conjunction with INT 
and/or WMAP. The values are quite robust as other experiments are added. (Similar 
values are obtained if we use the 1.8% Netterfield et al. (2002) cut rather than the 2.9%: 
0.953 ± 0.037, 0.938 ± 0.036, 0.949 ± 0.021 and 0.942 ± 0.020 for INT, INTc, WMAP and 
WMAP+INTc.) ACBAR was done with Boomerang included and the 2-14 bins. (The 
multipliers are lower when the first bin is included, 0.94 ± 0.05, 0.92 ± 0.05, 0.89 ± 0.047, 
0.914 ± 0.042.) ARCHEOPS was also done with Boomerang included for INT and INTc. 
The results are the same whether Boomerang is recalibrated or not. Beam adjustments 
are also computed by the method. For example, the Boomerang beam is compatible with 
no change, but substantially reduced errors: the prior is 1.00 ± 0.13 and we determine 
1.02 ± 0.03 with INT, 0.984 ± 0.026 with WMAP and 1.00 ± 0.022 with WMAP+INTc. 
When cosmological parameters are determined, adjustments such as these are done for 
each experiment to maximize the entropy of each model in turn, but the uncertainties are 
then marginalized. The optimal bandpowers also have the calibrations marginalized.) 
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ARCHEOPS 
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INT 
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WMAP+INTc 


1.00 + 0.10 
0.975 ± 0.035 
0.961 + 0.036 
0.963 + 0.017 
0.959 + 0.017 


1.00 + 0.1 
0.96 + 0.05 
0.947 ± 0.044 
0.912 ± 0.049 
0.938 ± 0.044 


1.00 ± 0.07 
1.06 + 0.034 
1.05 + 0.039 
1.046 ± 0.019 
1.046 ± 0.019 
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Table 2. Peak/dip locations £ pk /di P ,j an d heights C pk / dip j (in fiK 2 ) from maxi- 
mum-likelihood analysis of data within N bands as the bandwidth slides over the data. 
(Results for a four-band-width slider acting on the compressed grand unified spectra of 
the various epochs are shown. Values of the curvature \n p k/di P ,m\ in units of c Kpk are 
given beside the £ p k/di P ,j- A detection requires that this exceeds unity and £ p k/di P ,j lies 
within the bands probed. If a peak/dip detected with a three-band slider is more signif- 
icant than for the four-band slider, it is shown in brackets. When we do WMAP-only 
with this crude Boomerang binning we get the results in the second column that are in 
excellent accord with the Gaussian-fit results of Page et al. (2003): (220 ±0.8, 5373 ±337), 
(546±10, 2381±83) and (411.7±3.5, 1707±43) for the first and second peak and the first 
dip. Other individual experiments we have applied this procedure to include Boomerang, 
CBI and ARCHEOPS. The numbers in the table also accord well with those we have 
obtained for Boomerang (de Bernardis et al. 2002; Ruhl et al. 2003) and CBI (Pearson 
et al. 2003; Readhead et al. 2003) data alone; e.g. for Boomerang-only with the Faster 
power spectrum, recalibrated as in table 1, we get (221 ± 7, 3.5, 5371 ± 337), (528 ± 15, 
1.5, 2183 ± 136), (820 ± 15, 1.5, 2000 ± 215) for peaks and (412 ± 7, 5.2, 1706 ± 105), 
(699 ±21, 1.5, 1708 ±139) for dips. Benoit et al. (2003) report the ARCHEOPS first peak 
location to be at 220 ± 6; with our method, all sliders of width N > 4 detect the first peak 
accurately; e.g. N = 9 probing i values from 110 to 350 give (221 ± 6, 5.7, 5296 ± 360). 
The recalibration of ARCHEOPS of table 1 has been applied to C p k,i, bringing it into 
good accord with the WMAP and March 2003 numbers.) 
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Table 3. Peak/dip locations £ p k/di P ,j an d heights C p k/di P j (in f-K 2 ) determined by en- 
semble-averages over the Ce database and the weak prior. (Use of the weak prior allows 
large movement of peak locations associated with the geometry, hence is preferable to 
more restrictive priors for this application. These numbers should be contrasted with the 
'model-independent' numbers of table 2. The comoving sound speed, r s , and damping 
scale, Rd, and their associated angular scales, n£ a and Id, are also shown. Database num- 
bers using this method were given in de Bernardis et al. (2002), Pearson et al. (2003) 
and Ruhl et al. (2003) and the structural parameters for the June 2002 data were given 
in Sievers et al. (2003). The results here compare well with those given by Spergel et al. 
(2003) using the MCMC method: 220 ±0.9 and 535 ±2 for the first two peaks. The WMAP 
team also give r s — 144±4 Mpc for power-law ACDM models, and r s = 147±2 Mpc using 
all of the data and a running index. For -k£ s they get 299 ± 2 and 301 ± 1, respectively. 
When we do WMAP-only, we get similar values to March 2003 for the peak locations and 
amplitudes, but with larger errors beyond the third for the peaks that WMAP does not 
cover - errors more similar to those of the January 2003 data.) 
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Figure 4. ID likelihood curves for the weak+flat+LSS prior show that the parameters 
were in good agreement but also show how much WMAP has sharpened the picture. 
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Figure 5. Evolution of as determinations for the weak+fiat prior, with the tc prior applied 
as well for the March 2003 data. These are contrasted with the estimates from weak lensing 
and clusters which the LSS prior encompasses (Sievers et al. 2003; Bond et al. 2003b). 
The MCMC results give distributions that extend slightly more to higher as- 
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Table 4. Cosmological parameter values from the Ci database and their la errors are 
shown, determined after marginalizing over the other six cosmological and the various 
experimental parameters. (In all cases, the weak prior (0.45 < h < 0.9, age older than 10 
Gyr) is applied, and wq is fixed at —1, the cosmological constant case. The first set allow 
fl k to vary, the rest fix it at zero. The LSS prior agrees with weak lensing and redshift 
survey results, and agrees with most of the cluster determinations. The parameters are 
very stable if extra 'prior' probabilities for LSS are included, or if the HST range for h 
is used, or if SN1 data are included. Allowing wq to vary yields quite similar results. 
Although the optimal spectra include the January 2000 data in the subsequent mixes in 
figure 2, these parameter estimates do not. (Their inclusion has little effect.)) 
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Table 5. MCMC determinations of cosmic parameters should be compared with those of 
the Ce database given in table 4- (Means and standard deviations are given here, rather 
than the Bayesian 50%, with 16% and 84% 'la' errors. In addition, for the March 2003 
data, the Bayesian results using the 899 WMAP points in the database with the rc-prior 
are shown in the last column. Note how as increases when the logarithmic running of 
the spectral slope is allowed, a manifestation of the high correlation of as and rc with 
dn s /d\nk as well as n s . The values are therefore sensitive to the priors imposed.) 
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